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■ We analyze recent measurements of the room temperature current - voltage char- 

acteristics of gold nanowires, whose zero current conductance is quantized in units 
of 2e 2 /h. A faster than linear increase of current with voltage was observed at low 
voltages beginning from V c = 0.1 V. We analyze the nonlinear behavior in terms 
of a dynamic Coulomb blockade of conducting modes and show that it may be 
explained as a Luttinger-liquid effect. 
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1. Introduction 

Recently conductance quantization at room temperature was discovered jlj in metallic nanowires 
formed by breaking a contact between two metallic electrodes. Although the accuracy of the quan- 
tization in metallic nanowires as a rule is less then what can be achieved in ballistic nanocontacts 
formed in the two-dimensional (2D) electron gas of gated semiconductor heterostructures J2|, the 
mere fact that plateaus appear at quantized values of conductance G ~ n(2e 2 /h),(n = 1,2,...) 
shows that the electron transport through the nanowire is ballistic (the mean free path of charge 
carriers is larger than the length L and width W of the wire). In thin metallic nanowires the screen- 
ing of ionic potentials by mobile charge carriers, which is very efficient in the bulk electrodes, could 
be considerably suppressed. Hence, the electrostatic potential felt by electrons moving along the wire 
becomes more smooth than is the case in bulk metals. In this picture the conductance quantization 
at room temperature can be understood within the framework of the "conductance is transmission" 
- picture — originally proposed by Rolf Landauer |3| — applied to a ballistic contact with almost 
perfect transmission. The approach of Landauer was subsequently developed by him and others into 
what we now know as the Landauer - Biittiker formalism, which is what is used below (see e.g. the 
review in Ref. Q). 

Suppression of screening in metallic nanowires implies in addition that the effects of electron- 
electron interactions could be important for the nature of electron transport through a nanocontact. 
It is well known that for ID ballistic wires with a length that is effectively infinite, the interaction 
renormalizes the conductance quantum to a value which depends on the interaction strength (in the 



2 



Goteborg Preprint APR 97-22 



Tomonaga-Luttinger liquid model the conductance per spin orientation is K p (e 2 /h), where K p is 
the charge correlation parameter of the Luttinger liquid |^|, |^|). However, the conductance of a finite 
ballistic wire attached to leads with nonintcracting electrons takes on integer values (in terms of 
the conductance quantum e 2 /h) regardless of the interactions in the wire 0. Therefore, the effects 
of electron-electron interactions are not revealed in the linear response regime for a finite ballistic 
(reflectionless) contact. 

If the electron transport through a wire is not purely ballistic (i.e. there is a finite probability for 
electrons to be backscattered by impurity- or channel imperfections) the above picture is radically 
changed. In an effectively infinite ID wire even arbitrarily small backscattering suffices to block the 
dc-current at vanishingly low voltages and temperatures ^ [| . The effect is due to a strong enhance- 
ment of the bare backscattering amplitude by charge fluctuations at long distances. For a finite wire 
the renormalized barrier is also finite; the procedure of barrier renormalization should be stopped 
when a low-energy scale of order e 2 / L is reached. Therefore, in a multichannel quantum wire the 
transmittivity of quantized modes will depend both on their bare transmission coefficients and on 
the strength of electron-electron interaction. The transmitting modes most sensitive to backscatter- 
ing (the criteria will be given below) are converted to reflected modes by interaction effects. In other 
words the renormalized barrier is the actual barrier experienced by electrons propagating along the 
wire and this very barrier separates transmitting- from reflected quantized modes. The blockade 
imposed by the electron-electron interaction is gradually lifted by the bias voltage beginning from 
V > V c — e/L resulting in a power-law I — V characteristics. 

Nonlinear I — V characteristica for metallic nanowires in a regime where the linear conductance 
was quantized have recently been measured by Costa-Kramer et al. JTo[ . The following remarkable 
universal properties were observed: 

(i) the characteristic crossover voltage V c ~ 0.1 V to the nonlinear regime is sample-independent 
and anomalously small in comparison with the atomic energy scale ~ 1 eV set by the spacing 
of quantized mode energies in the contact pd[ 

(ii) the nonlinear character of the I — V curves does not depend on the conductance in the zero 
voltage limit and is therefore independent of the number of conduction channels 

(iii) the nonlinear contribution to the current is approximately proportional to the third power of 
the voltage 

The purpose of this work is to show that all the above listed properties can be explained (at least 
qualitatively) if one associates the nonlinear current with the current produced by one (or possibly 
a few) — partly reflected — modes in a multichannel wire of interacting electrons. The current 
associated with this mode, which is small due to "Coulomb blockade" effects in the linear response 
regime, is strongly enhanced when the bias voltage is increased; eventually the voltage lifts the 
blockade and converts the reflected mode to an almost perfectly transmitting one. As we proceeed 
to elaborate this conjecture we recall first how the universal power-law dependence of current on 
voltage arises in the Luttinger liquid approach. To be concrete we consider a ID model of interacting 
spinless electrons scattered by a local potential barrier. 

2. Luttinger liquid model for nanowires — Weak interaction limit 

For our purpose it is important to be able to describe a situation where the bare transmission 
probability can be close to — but not necessarily equal to — unity so that the conductance increases 
by essentially one quantum when the Luttinger liquid effects are suppressed (for example by a high 
bias voltage) . This opening up of a new conduction channel is precisely what we propose may explain 
the nonlinearities in the measured I — V curves. For arbitrary barriers the model can be treated 
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analytically only for weakly interacting electrons when the effect of interaction can be reduced to 
a barrier renormalization produced by the formation of Friedel oscillations of the electron density 
near the barrier. The renormalized transmission coefficient is fl2| 
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Here To(i?o) is the bare transmission (reflection) coefficient (To + Ro = 1), £f is the Fermi energy, 
-D is an ultraviolet cutoff (see below), 7 = 7/ — 7b = [V ee (0) — V ee {2kF)]/2irTiVF is a dimensionless 
parameter which characterizes the strength of the electron-electron interaction, and V ee (q) is the 
Fourier transform of the interaction potential. Strictly speaking the parameter 7 should be small for 
this approach to be valid (7 ~ e 2 /2itTivf <C 1), but it is believed that Eq. (f ) still holds qualitatively 
in the range 7 ~ f . In what follows we will regard 7 as an input parameter. 

The renormalized transmission coefficient (1) leads to an expression for the current when put into 
Landauer's formula. At zero temperature one finds 
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We assume that the bare transmission amplitude, To, is energy independent near the Fermi level, 
so that the entire energy dependence of T R (e) stems from the universal effects of barrier renormal- 
ization. The function F 1 {x) has to be obtained numerically except in special cases (see below). 
For a small renormalized transmission probability (T R <C 1) one can approximate Eq. (1) as 
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The physical meaning of the high-energy cutoff Do in Eqs. (1) and (3) is the maximum energy 
transfered in a scattering event. For a strictly ID wire with a local electron-electron interaction 
(Luttinger liquid) a natural cutoff is the bandwidth £p- If the interaction potential is characterized 
by a finite length scale l s (a screening length) one has instead D ~ Hvf/I s ■ In a thin (l s ^,W) 
nanowire screening effects could be suppressed. In this case the Coulomb potential averaged over 
quantized transverse modes takes the form V c (x) ~ e 2 j\/x 2 + W 2 and the relevant high-energy 
cutoff is reduced to Do ~ Tivf/W . The physics at scales x <§; W can not be described by a long 
wavelength approximation and the Luttinger liquid-like approach is inappropriate. 

The reasoning outlined above leads to an expression for the nonlinear tunneling current which 
takes the following form (eV <§C Do) 
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We will model a nanocontact formed in conductance quantization experiments as an atomic-size 
constriction (characterized by a set of transmission coefficients for quantized modes) attached to a 
long (L 3> W) nanowire of width W (see Fig. I). In an adiabatic model of a nanoconstriction [[l3| 
it is reasonable to assume that To is exponentially close to zero or unity for all channels but one (or 
possibly a few) . It should be stressed that in the range 7 ~ 1 Eq. (4) holds even for To close to unity 
due to a strong renormalization of the transmission amplitude. The crossover value for the ratio 
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Fig. 1. Model of a nanowire contact. An atomic size constriction (within dashed circle) characterized by a set of 
transmission coefficients for quantized modes is attached to a long nanowire of of length L and width W . Most modes 
are either fully transmitted (lower arrow) or completely reflected (not shown) by the constriction. One, or possibly a 
few modes are partially reflected (upper arrow). 



of unrenormalized coefficients Tq/Rq that separates "under-the-barrier" from "over-the-barrier" 
transport is {Tq/Rq) c — (L/W) 2 ' 1 which is much bigger than unity if 7 ~ 1. If To is close to one we 
should use the exact expression (I) rather than Eq. (3) when describing high voltage behavior. 

Since the nonlinear contribution to the currents measured in Ref. |1(J was fitted to a power law, 
V q , with an exponent q statistically peaked at q ~ 3, the specific value of interaction strength 7 = 1 
appears to be of most interest to us. In this special case the nonlinear conductance calculated from 
Eqs. (1) and (2) takes a simple analytic form 



Therefore, if the "Coulomb blockade" of the quenched tunneling mode is lifted, the nonlinear voltage 
dependence of the current goes over to the linear dependence expected for a perfectly transmitting 
mode. As for the modes which are propagating, "over-the-barrier" reflection only results in a small 
additive renormalization of the linear current. 

Is it reasonable to trust the power law I — V result even in cases when the coupling strength is 
not small? We think so; actually this is a quite general result for Luttinger liquid-like theories (see 
e.g. the reviews H]). The analogous expression to Eq. (4) can be derived in the tunnel Hamiltonian 
model for the whole range of interaction strengths. However, in this approach one has to assume 
a small unrenormalized transmission coefficient (To < 1) to get analytic results |L For our 
purpose it was important to describe scattering of interacting electrons by a barrier of arbitrary 
height, therefore we have used the approach elaborated in Ref. p~2| . 

For a finite wire of length L renormalization of the transmission amplitude should be stopped 
at energies Se = \e — Ef\ _ A = litfivplL . In the weak coupling limit the low-energy cutoff A 
corresponds to the minimum-energy charge excitations in a finite wire (longitudinal energy dis- 
cretization). In the case of strong interaction the cutoff depends on the plasmon velocity s and is 
A(s/vf) 2 ~ e 2 /L (see Section 3). When describing the current in a finite wire Eqs. (4) and (5) 
are valid for voltages V > V c ~ A/e .At lower voltages the transmission coefficient is energy 
independent T R (8e ~ A) = (To/ Rq)(W/ L) 21 . The current produced by the quenched mode at low 
voltages, V < V c , is linear and small even if the bare transmission amplitude is close to one, 
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2.1. Comparison between theory and experiment 

Employing the above simple formulae for an explanation of the nonlinear conductance of metallic 
nanowires we first note that the height of the measured conductance steps deviate somewhat 
from multiples of e 2 /h (i.e. from perfect quantization). By associating the difference 5G with the 
small current carried by Coulomb blockaded mode at low voltages we conclude from Eq. (6) that 



5G/(e 2 /h)~(T /R )(W/L) 2 ~<. (7) 

As discussed above the suppression of the current in this quenched mode is lifted by a sufficiently 
large applied voltage. This gives rise to a nonlinear contribution to the I —V curves when V ~ V c . It 
follows that the nonlinearity is due to a suppression of Luttinger liquid-like effects at high voltages. 
Since the nonlinear contribution to the current arises from the mode most sensitive to backscattering 
it does not depend on the number of conduction channels (see below). The universal character of the 
nonlinear contribution is due to the "long distance" origin of barrier renormalization; the relevant 
energy scale is either given by the total length of the nanowire L (if V < Vc ) or by the "field" 
length l E = hv F /eV > W. 

Simple numerical estimations show that the above formulae allow us to explain the I — V char- 
acteristica measured in Ref. First of all the crossover voltage observed in the experiment 

(V c — 0.1 V) can be readily obtained from Eq. (5) for metallic nanowires a few nanometers long. 
Since the numerical coefficient in the theoretical estimation of V c is uncertain it is difficult to be 
more precise in the quantitative predictions. However, it is worthwhile to stress that the voltage in 
question is determined by the total length of the nanowire and not by the length of its narrowest part 
which provides the conductance quantization. Therefore the crossover to the nonlinear regime in 
this picture does not depend on uncontrollable factors like precisely how the constriction is formed. 

From the experiments of Ref. Jl0| we conclude that the exponent (2j + 1) in Eq. (4) should be 
approximately equal to 3, hence 7 ~ 1. For Coulomb interaction V ee (0) 3> V r ee (2fcF) and one can 
estimate 7 ~ (e 2 /itTivf) \n(L/W) which for reasonable values of L/W and vf typical for metals is 
in amazingly good agreement with the value extracted from experiment. This numerical coincidence 
could be a strong argument in favour of our ideas if both theoretical predictions and the chosen 
"experimental" value of exponent (q — 3) in power-law I — V curves were quantitatively reliable. 
Unfortunately there is no physical reason to believe that Eq. (1) is quantitatively correct in the 
range 7 ~ 1. As for experiment [fiof , the distribution of the exponent q for different measurements, 
although peaked at q = 3, spreads out from 2 to 4. So instead of playing with numbers it would be 
more reasonable to look for qualitative features which can be easily either confirmed or ruled out by 
experiments. 

The qualitative predictions which emerge from our picture are straightforward to extract from 
the theory. First, the larger the fluctuations of the conductance around its quantized value the 
bigger the nonlinear current. One can find some experimental confirmation of this prediction by 
examining Fig. 11 of Ref. jn^. (The fluctuations of conductance are obviously increasing as the 
number of conduction channels [modes] increases. So one can expect that the absolute value of the 
nonlinear current will increase as well. Such a tendency can indeed be observed in Figs, ll(a-d) of 
the cited paper). Next, it is evident that the contribution of the tunneling mode to the current when 
the Coulomb blockade is lifted by high voltages could at most be comparable to that of a freely 




Fig. 2. Conductance G(V) = I /V calculated from the Tomonaga-Luttinger liquid model, Eqs. (5) and (6). It 
represents the contribution of a conductance channel in which electrons are partly reflected (To < 1) from a potential 
barrier renormalized by electron-electron interactions, the strength of which is characterized by a parameter 7. The 
parameter Dq = hvp/W is a high energy cutoff associated with the time it takes for a charge excitation to spread 
across the width W of the wire, while V c ~ hvp/L is a low energy cutoff related to the finite wire length L. 



propagating mode (for which Rq -C 3b). Therefore the total current of an n-channel constriction at 
saturation voltage (which is determined by a short distance length scale) should approach the linear 
current of an (n + 1) channel nanocontact. 

Plots of the nonlinear conductance G(V) calculated from Eq. (5) are shown in Fig. 2. Even for 
a very nearly ballistic mode with To = 0.9 one can see that there is a significant (although not 
complete) restoration of transmittivity for voltages in a range relevant for the experiment. However 
for the quenched modes with small backscattering complete restoration is achieved at high voltages 
V ~ V s ~ y/To/ RqDq ^> Dq, where the Luttinger liquid approach does not hold. In the theoretically 
reliable region of voltages eV < D the current carried by each individual (quenched) mode is a 
few times smaller than for a perfectly transmitting channel. In the experiments [fl0|| , on the other 
hand, the conductance is increased by approximately one quantum unit already in the low voltage 
region. In our picture the observed conductance enhancement could be explained if more than one 
suppressed mode contribute to the nonlinear current. 

2.2. Spin degrees of freedom and multichannel effects 

Although the above formulae pertain to the case of spinless fermions, the inclusion of spin does not 
change the situation radically. As was shown in Ref. Jl2| , adding the spin degree of freedom only 
results in an obvious additional factor of 2 in the conductance for the case of a smooth interaction 
potential, V ee (2kp) <C V ee (0). In a general case the transmission probability in Eq. (1) should be 
replaced by jl2l 
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It is evident that the additional logarithmic energy dependence of the transmission probability can 
not change our estimations in an essential way. 

A less evident point in connection with the above purely ID considerations is how they are affected 
by multichannel effects. In reality a metallic nanowire may only be able to support a few conductance 
channels at its most narrow part where the width is of the order of a few Fermi wavelengths. The wire 
gets much broader, however, away from this region and many conduction channels are supported. 
To discuss this complication we will first argue that even though the potential barriers that define 
the wire are due to atoms the adiabatic approximation applies. 

In principle a rearrangement of the atomic structure as the the wire is extended could lead to a 
step-like behavior of the conductance, but there is no reason to expect such conductance steps to 
be quantized in units of e 2 /h. Experiments, on the other hand, have definitely demonstrated such a 
conductance quantization (statistically) of metallic nanowires (see e.g. the review Jig]) . Moreover, 
the conductance staircase measured by the STM technique at room temperature |10| looks very 
similar to the one found in 2D electron gas systems. In particular, the accuracy of the quantization 
is quite high (within 10% for the first step) and all plateaus were observed (up to n = 6 — 8) with no 
indication of suppressed conductance steps at n = 2,4, 7 as predicted for a cylindrically symmetric 
ballistic microcontact. This means that for the selected measurements, and when the conductance 
is quantized, the motion of electrons through the constriction can be considered adiabatic; so a 
simple "2D" -model (with no degenerate transverse modes) of an adiabatic contact jl3| may be quite 
adequate for our purposes. 

In the adiabatic approach to electron transport through a quantum contact the transverse (fast) 
and longitudinal (slow) motion of electrons through the contact are decoupled when the Schrodinger 
equation is solved. The kinetic energy due to the transverse motion then appears as a potential bar- 
rier for the longitudinal motion which adds to the smooth potential from the constriction itself. 
The smooth effective potential is different for different modes and has a maximum where the con- 
tact diameter has its minimum. Since the energy level spacing imposed by transverse momentum 
(mode) quantization is of the order of the Fermi energy, one can use the following simple classi- 
fication of modes: those with energies above the barrier are transmitting modes, those below the 
barrier reflected modes. As we have seen already this classification should be changed for interact- 
ing electrons, where the renormalized potential rather than the bare potential is relevant. One can 
infer from Eq. (1) that all modes for which Tq/Ro -C (L/W) 21 ^> 1 are quenched by the strong 
enhancement of backscattering. It is known that for smooth potentials over-the-barrier reflection 
decreases exponentially with increasing energy of the incident particle. Hence it is quite natural to 
assume that only the last mode (the one with energy closest to the potential barrier maximum) will 
change its character by renormalization effects and be converted from a transmitting to a reflected 
mode. All other modes keep their unrenormalized character; they are either transmitting modes 
(over-the-barrier) or totally reflected modes (under-the-barrier). 

Let us now ask the question how a transmitting mode with T = 1 is affected by inter-modc 
(inter-channel) forward scattering in a multichannel wire. The answer is that in an effectively infinite 
multichannel Luttinger-liquid wire the Coulomb interaction smears out the quantized conductance 
plateaus so that in the strong interaction limit the conductance is described by almost smooth 
function of the wire width (or of the Fermi energy) |l8| . This effect is entirely due to the properties 
of plasmon excitations at long distances. For a wire of finite length connected to large reservoirs 
the dc current will be determined by the properties of plasmons in the reservoirs [p~8L [7fl. So in 
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experiments, where the reservoirs are massive metals, the conductance of a perfect multichannel wire 
of interacting electrons will be quantized exactly as if the electrons in the wire were noninteracting. 

Now we turn to another question: Could the Coulomb-blocked mode, which is reflected by the 
renormalized barrier, affect the propagation of the transmitting modes? In the case of weak inter- 
action the answer is that it definitely cannot. First of all the presumed smooth bare potential can 
not induce intermode scattering. Therefore each mode with ]?J ^ (j is the mode index) sets 
up its own Friedel oscillations when reflected by the potential. The corresponding induced barriers 
are characterized by different oscillation periods \(s — e F ) , where e F is the "Fermi energy" 
in the j:th channel. (The channel dependence comes from the fact that e F is measured from the 
quantized transverse energy corresponding to the j:th mode; it is the bandwidth of the j:th mode). 
Interaction enhanced backscattering arises only for electrons in the same mode or subband [since 
the corrections to bare transmission amplitude is proportional to Rq ]■ Particles with energies in 
the vicinity of e J F are phase locked to the Friedel oscillations and thus strongly backscattered. 
The influence of this effect on the transmittivity of quantized modes have already been studied 
above. Perfectly transmitting modes (Rq^ — 0) can be scattered by the potential due to Friedel 
oscillations of the charge density in the quenched channel. However, this scattering could lead at 
most to nonsingular (length independent) contributions. In the lowest order of perturbation theory 
the corresponding induced backscattering amplitude \r^' \ — 7r-y | ro | ( ro is the bare backscattering 
amplitude of quenched mode) is small if \r \ <C 1. 

To put it differently, transmitting modes will propagate freely even in an interacting electron 
system irrespective of possibly strong backscattering in the other modes. Is it possible to extend 
this claim to the strongly interacting case? To clarify this point we switch now to the case of strong 
coupling, where the picture of Luttinger liquid transport through a constriction will look slightly 
different. 



In the strong coupling regime the suppression of tunneling for repulsively interacting electrons can 
be easily understood in terms of a Wigner crystal (or charge density wave) pinned by impurities Q . 
The quantum depinning of a Wigner crystal is accompanied by propagation of plasmon excitations 
along the wire. The gapless sound-like spectrum of ID plasmons results in a logarithmically divergent 
contribution of plasmons to the tunnel action. The infrared cutoff (at small momenta) appears at 
scales determined by temperature or voltage. However, in the case of a finite wire (either isolated 
or connected to 3D leads) one gets the additional natural energy scale A s = Tis/L, where s is the 
plasmon velocity. For an isolated wire (or a ring, see e.g. Jl7]]) this scale is provided by the discrete 
spectrum of energy levels; for a wire connected to 3D leads a cutoff appears due to the fact that the 
excitation spectrum of 3D plasmons has a gap and does not give rise to infrared divergencies. 

As in the case of weak coupling we will associate the nonlinear contribution to the current with 
the interaction-quenched mode. If the interaction is not weak one can treat the Luttinger liquid 
transport through a barrier analytically only in the limit of a strong barrier. So our consideration 
could be applied strictly speaking to the tunneling regime of conductivity. We start with a purely 
ID picture (a single mode channel). Now the current at temperatures effectively equal to zero (for 
metals this temperature interval includes room temperature) is given by 



3. The strong coupling regime 




(9) 



where 
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is the tunneling probability associated with the excitation of long wavelength plasmons ]^] (notice 
that in this reference the exponent in the formula for the tunneling probability contains an incorrect 
numerical factor). In Eq. (10) V p is the pinning potential which for strong repulsion (a <C 1) can be 
put equal to the height of the potential barrier ||. 

At high voltages eV > A s /a (however still much lower than the height of the barrier) the current 
is nonlinear 

I SL (Y)^aV p (^y . (11) 

Notice that we omitted terms of order unity in the exponents of Eqs. (10) and (11) since the accuracy 
of calculations in the strong coupling limit (a <C 1) is not so good as to make it meaningful to keep 
them. 

At small voltages the tunnel action should be cut off at the energy scale A s and the tunnel current 
becomes linear 

To get the kind of nonlinear I — V characteristics measured in Ref. [H one should put a ~ 2/3 
in Eq. (12). Again, this value of the correlation parameter is not sufficiently small to justify a 
quantitative comparison between theory and experiment. 

An equation for the nonlinear current [analogous to Eq. (11)] valid in the required range of 
interaction strength can be derived in the tunnel Hamiltonian approach fig] . The corresponding 
formula for a finite wire connected to ID reservoirs of noninteracting electrons (for ballistic transport 
an analogous model was considered in Ref. 0) takes the form 



4e \t\ 2 f°° ( eV 2 / A \ 

/(V) = — V- / dxcos< — x arctan — x ) + (a^ 1 - l) \tt + 2si(x)} 

h A s J \A S a \A S J y n 

exp {2 (a" 1 - l) [C + ln(x) - ci{x)] } , (13) 

where t is the tunneling matrix element which is assumed to be small, C ~ 0.577 ... is the Euler 
constant, A is the ultraviolet cutoff (A > A s , A ^> eV), and the integral sine (cosine) function is 
defined as follows 



1 ' h 



dt 



sin 

cos 
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(14) 



For A s — * (infinite wire) Eq. ( J13|) yields the well-known result for Luttinger liquid transport of 
a power-law I — V dependence P,]l5[| 
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Fig. 3. Log-log plots of the nonlinear tunneling current (Eq. 13) vs. gate voltage for a finite quantum wire 
connected to reservoirs of non-interacting electrons. Here Iq = (e/h)jt| 2 /A, where t is the tunneling matrix element, 
A ~ hvp/W is a high-energy cut-off related to the wire width W; A s = hs/L (s is the plasmon velocity) is an energy 
scale associated with the discrete energy spectrum and hence the finite wire length L. In the calculation A/A s = 50. 
For small voltages, V < A,/e, the current increases linearly with voltage. For higher voltages the current depends on 
the interaction parameter a and the slope in the figure approaches the infinite-wire length value 2/a — 1 of Eq. (15). 



e \t\ 2 feV 



1JVL — T a i i 

h A V A 



(15) 



For a finite wire the integral in Eq. ( |l3| ) has to be evaluated numerically. The corresponding I — V 
curves are plotted in Fig. 3. From this figure one can see that at low voltages the current grows 
linearly with voltage and is small. Beginning from V ~ V c \ ~ A s /e the tunnel current increases 
rapidly, eventually reaching the asymptotic value Eq. ( |l5| ) at V > V c i ~ A s /ea. 

As we have noticed already the above formulae describe the tunnel current and strictly speaking 
they can not be used for explanaining the nonlinear I — V behavior measured in the experiment 
[^p| (where the observed increase of conductance was of the order of the conductance quantum). 
However, in order to analyze the phenomenon in terms relevant for the strong coupling regime 
one may put \t\ ~ A ~ hvp/W . For a comparison with experiment it is reasonable to use the 
expression [analogous to Eq. (jl5|)1, which describes the nonlinear current of spin-1/2 electrons. The 
corresponding formula can be readily obtained from Eq. (|l5| ) by the substitution loT 1 — » 1 +a o r 1 , 
where a p (a a ) is the Luttinger liquid correlation parameter for the charge (spin) sector). For an 
SU (2) invariant spin interaction (relevant here) a a = 1 . The correlation parameter in the charge 
sector can be expressed in terms of the Fourier transform of the forward scattering interaction 
potential (see e.g. Ref. ||i"8f ) 



(16) 
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For a cubic I— V behavior a p ~ 1/3 and the corresponding value of the dimensionless interaction 
strength is 7 ~ 2 , which is larger by a factor 2 than the value extracted in the weak coupling 
limit. All numerical estimations performed for the weak coupling limit can be equally well reproduced 
using formulas corresponding to the strong coupling regime (obviously with the same uncertainties in 
numerical coefficients of order one). As a result one has to admit the evident failure to quantitatively 
describe nonlinear transport in metallic nanowires using naively a Luttinger liquid model. However 
it is suggestive that both limits yield reasonable values for the interaction parameters (of order one, 
which is typical for metals). 

The advantage of looking at our problem in the strong coupling limit is that it allows us to consider 
the multichannel case analytically. At first we assume that the current through the nanoconstriction 
is due to tunneling. Then the problem is reduced to describing the propagation of charged excitations 
(plasmons) in a multichannel wire. For this case (here we again for simplicity will consider spinlcss 
fermions) it is known Jl8| that in the limit of a large number of channels N± ^> 1 the plasmon 
velocity s is given by 



\ 



N ± 2 

: f m 
2=1 



where nj is the electron density in the j-th channel. Hence it is determined by the total density 
of electrons in a multichannel wire. It is this very density (which approaches the bulk density 
of electrons in the vicinity of the electrodes) and not the smaller density of electrons in the few 
quantized modes at the center of the contact which will determine the actual value of the correlation 
parameter a in Eq. (11). Therefore, the exponent in the power-law type I — ^-curves should be 
independent of the number of conductions channels as was claimed in Ref . jlO| . 

A situation where the number of channels varies along the wire can be modelled by assuming 
the plasmon velocity to be a function of position s — ► s(x) = s f(x). Straightforward calculations 
analogous to those performed in Ref. || (see also |^7| ) yield 

£(£) = exp(-— dx^ \ . (18) 

\ "0 Je x J 

According to Eqs. (9) and (18) the nonlinear conductance of a general multichannel wire will depend 
on the actual shape of the wire. However, if the function s(x) is smooth, the tunneling probability 
still has a local energy dependence, 



D(e) = (j^ ° <e> , a(e) = v F /s{e) . (19) 

Therefore to get in this approach a sample-independent exponent describing the nonlinear conduc- 
tance observed in Ref. JlCj ] one has to assume that the nanowire has an almost uniform shape far 
from the constriction (this is exactly the model for a nanocontact used here, see Fig. 1). 

The last point we would like to discuss is how the mode quenched by interaction effects influences 
the transmitting modes. In the strong coupling regime the assumption of perfect transmission means 
that the corresponding Wigner crystals can slide freely (without distortions) along the channels. If 
the interaction potential is smooth and one can neglect effects of backscattering [which are propor- 
tional to the exponentially small interaction constant V(2kp)] the degrees of freedom associated 
with perfectly transmitting modes decouple from the rest. Their dynamics is linear and can be 
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treated exactly (7| for an adiabatic constriction. One may conclude that as long as the backward 
scattering part of the interaction potential is irrelevant, perfect conductance quantization can not 
be spoilt by quenching of reflected modes. 



4. Conclusions 

Summing up all facts in favour of the above picture, we conclude that the Luttinger liquid-like effects 
could be a reasonable explanation of the nonlinear conductance measured in metallic nanowires. If 
so the experiment by Costa-Kramer et al. |Io|] can be regarded as one of only a few highly important 
experiments (see also |H| |2Cj ) where the fundamental properties (Kane-Fisher effect) of Luttinger 
liquid-like states of strongly correlated electrons were observed. 
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